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1.  INTOUllllCTlON 


Given  n set;  9C  {x  : Ax  ,<  a)  and  a J.lnear 

function  /{(y)  ■ Cx  + c from  to  r”  wo  consider  tho  cKlstenco 
and  computation!  In  a finite  nunijor  of  atepa,  of  a atatlonary  polnti 
A point  X*  in  la  defined  to  be  n stationary  point  of  ((fJK)  if 
Cx*  + c la  an  inward  normal  of  X Ht  x*,  or  in  other  wordu*  If 
(y  •<  X*)  • ^(x*)  la  nonnogatlve  for  all  choices  of  y in  9^  • 

The  existence  and  computation  of  stationary  polntu  is,  in 
particular,  central  to  tho  solution  of  certain  quadratic  prn&rama, 
matrix  gamoa,  and  economic  cquiUbrluin  problcinsi  Such  problems  can 
often  be  caat  into  the  linear  complementarity  problem  which  is  a 
atatlonarlty  problem  (^,3f)  where  Is  tho  nonncgatlve  orthant; 

Lemks's  algorithm  [6]  offers  tho  principal  avenuo  for  solving  the 
linear  complementarity  problem. 

Herein  we  adapt  Lomko's  algorithm  in  order  to  approach  tho 
general  statlonarlty  problem  Towards  deacrlblng  our  main  re- 

ault  let  x^  be  an  arbitrary  point  in  X • Wo  introduce  constraints 
Bx  .<  b so  tliat  x^  Is  the  unique  solution  to  Ax  < a and  Dx  < b. 

Next  define  to  be  the  aet  of  x*s  such  that  Ax  < a ond 

Bx  £ b + Oe  whoro  e ■ (1,  1,  1)  and  0 varies  over  the  in- 

terval [0,  ■f»)  • Alao  define  a piocowlae  linear  path  (X,0)  to  bo 
a (continuous)  function  (X,0)  t [0,  +»)  R**  x tO,  +»)  that  is  w 

Whi(»  Ststlen  If 

offine  on  sacli  of  a finite  number  of  closed^lntervalu  where  the  closed  Buft  SittlM  g 
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Intervnlii  covor  [t),  4'“)  . The  following  thooi  em  iMiMureii  our 
pritu'lpnl  conclusion. 

Thcoroin!  Tlio  nlgotithm  computes  a piccBwlae  linear  path 
(X,®)  such  that  X(0)  - x®,  0(0)  ■ 0,  0(p)  tends  to  Infinity  as 
p does,  and  X(p)  Iti  n stationary  point  o f for  nil  p.M 

Ono  interest! nr,  consoquonco  of  the  tlieorcm  is  the  following. 

* f*'®  alsorlthm  computes  either  a stationary 

point  x''  of  (t,X)  ^ >‘‘'y  + Ox  j 0 ^ 0*}  In  3C  such  that 

X • + Ox)  is  negative  for  all  0 exceeding  O*.  B 

for  emphasis  and  clarity  wo  note  thot  throughout  the  paper 
C is  not  assumed  to  be  aymetrir.t  K Is  hot  aosuinod  to  lie  in  the 
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In  this  section  v/o  give  a dascrlptlon  of  tho  iilgorlthm. 

Certain  details  of  exocutlon  arc  ampllflod  upon  In  Section  3,  a 

t 

complete  example  is  executed  in  Section  4,  and  the  proof  of  con- 
versance Is  elaborated  upon  in  Appendices  1 and  2. 

Using  a theorem  of  the  altornatlve  it  Is  A slmpla  matter  to 
eea  that  finding  a stationary  point  x of  is  equivalent  to 

finding  a solution  to  tho  system 

T 

Ax  + 8 “ ii,  Cx  + c + A X • 0» 

(1) 

If  X is  a stationary  polnti  then  tliere  is  a (s»  X)  no  that 
(x,  S|  X)  la  a solution,  If  (x,  s,  X)  ia  a solution  then  x la 
a stationary  point. 

T 

With  regard  to  the  data»  A is  an  m x n matrix,  A la  the 
transpnae  of  A,  a Is  m x 1,  C ia  n x n,  and  c is  ns],  and 
for  the  variables,  x la  n x i,  s Is  m x i,  and  X is  m x i. 

Step  One  of  the  algorithm  Is  to  select  any  point  x®  Inflf.  H 
In  Sactlon  3 ve  describe  a way  of  executing  Slop  One,  If  It 
Is  noted  that  % it  a aingletoni  then  the  algorithm  could  be  terminated 
for  x^  is  e stationary  point i however,  it  is  not  necessary  to  terminate 

I 

here  in  this  case,  . 
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,IV'’  aljiorl  t liiM  Im  to  atjjoln  mhl  1 r loiial  <'OniU  r.i ini  s 

Bx  < b In  ofdor  that  x^  ifl  tlio  »ml<|iui  Holutloii  to  the  nyHicm, 

Ax  < a,  Bx  j<  b • ®i® 

- * 

Aesutn:lng  I additional  conatralnta  are  added,  that  la,  (B,  b) 

in  it  X (n  -f  1),  tlicn  obncrvo  that  I muat  bo  at  leant  an  ltirp,o  an, 

but  need  not  bo  larpor  than,  n - li  + 1 whora  h In  the  nutiibor  of 

llnonrl/  Indepenflont  rowa  of  A nurh  tlint  A^^  . By  ( 

wo  denote  the  1th  row  of  the  oncloned  matrix.  In  Section  3 wo  dencrlbc 
a way  of  nuluctlng  (B,  b). 

Towordn  flmling  a atatlonory  point  of  wo  conulder  the 

Intermodlnry  taaU  of  finding  atatlonary  polntn  of  (t*  J^Tq)  where  JCq 
la  the  aet  of  x'a  auch  that  Ax  ^ a and  Bx  ,<  b ^ Oe  where 

e “ (1,  ....  1).  Clearly  Xq  empty  for  negative  0 and  boundod 

for  all  0. 

Computing  a atationary  point  of  (g,  3Cg)  ia  equivalent  to 
solving  the  system. 

Ax  + s"a,  Cx  + c+  A^X  + b'^'p  “ 0, 

(2)  Bx  + t - 00  ■ b*  b*A  ■ t'M  “ 0, 

■iO,  t>0,  ^iO,  hiO. 

Of  course,  if  we  have  a aolucion  (x,  a,  t,  X,  p,  6)  to  (2)  with 
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p ■ 0 tlicin  (X|  B,  X)  iiolvcB  (1)  nnil  x is  a stationary  point  of 

(J.». 

RawrltJns  (2)  in  detachod  cooriiicliMiL  Dorm  wq  bavo 


X s t X M''-  . • 0 


[Z 

0 

n 

■1 

m 

0 



n 

I 

0 

0 

0 

OI 

m 

0 

1 

0 

0 

umaiii 

Dll 

a > 0|  t i 0,  X >.  0,  P i 0,  0 _>  0, 

# • X - t • p * 0 . 


For  purposQH  of  computational  effic fancy  and  to  emphaala  the 
relation  to  I.cmke's  olgorithmt  we  now  oliminuto  x from  the  system 
(2).  By  ( ),g  we  denote  the  submatrix  of  columns  indexed  by  0i 

Step  Throo  of  the  algorithm  ta  to  select  a sot  8 of  n indices 
J from  fl.  m + 1}  so  that 


j e 3 


so  that  the  rows 


are  linearly  Independent, and  so  that 


(^,li)p  0 


hfl8  a oolution. 

In  Section  3 we  describe  a wny  of  aelectlnf,  0> 
Poi'  any  x In  JCo  w®  hnVo 


Stu|)  ,Ko«,r  of  the  ulj'.orlthm  1h  to  use  tho  expression  (3)  to 
eliminate  x from  tho  syiuom  (2)  In  order  to  obtain  the  system  (4) 
which  we  have  Illustrated  in  detached  coofflcient  form. 


H 

l 

X 

V 

0 

‘^1 

'^2 

A^' 

‘>1 

'^3 

0 

0 

‘>2 

s > 

0,  t > 

0,  X > 

0»  n > 

. 0,  s • 

• U ■ 0 . 


IiBttlnQ  a be  those  elements  in  {1,  ti  4-  1}  not  in  0 


we  have 
J 


etc.  {3 


Honco  (x,  9t  t,  A,  u)  noivofi  (2),  If  mid  only  If  Jt  iin.l.vpH  (:)) 
and  (4).  Thus,  if  wo  hoIvc  (4)  nnd  uno  (3)  to  compute  x,  we  Imvi* 
aolvod  (2)i 

Aa  la  done  in  the  simplex  metliod  In  order  to  remove  coincidence 
and  consequential  amblgulcyt  we  perturb  the  vectors  and  1^'’ 

qj^  + (a\  fi”'*  ‘I2  ^*'■^2 

2 n 

where  t la  a positive  infinitesimal,  [e]j  “ (e,  e , c ) nnd 

(c)^  ■ The  system  (4)  la  thus  porturbtid  to  ob- 

tain (5). 


(5) 


a 

t 

X 

U 

e 

‘^l 

a'*’ 

T 

“1 

'3 

0 

0 

**2 

8^0 

» t i 0 

• A ^ 0 

1 U 1 0, 

® i 0 

a • A ■ t • u “ 0 . 


That  la  to  say,  (3)  and  (5)  are  equivalent  to  (6). 
Equivolenuly,  defining  [ej^  and  (e)^  by 


<2)  la  perturbed  to  (6) 


7 


X 


H 


t 


X 


0 


coo  a'^’  0 -cTJa’S  In], 

A 1 0 0 0 0 fl+felj' 

no:00-e  b+  tc]^ 

B > 0,  t .>  0,  X >.  0,  u 1 0,  0 >.  0,  a • X - t • n - 0 . 


We  observe  thot  (a^i  t®,  » U®  i 0)  la  a solution  to  (5) 

C*  b C V 


where 


*0 


.y.  V7, 


+ c)  + [t], 


(■:•  •:), 


(■:•  -5. 


This  solution  of  (5)  Is  defined  to  bo  Iho  initial  solution.  We  shall 
show  in  Appendix  1 that,  in  fact,  tho  Initial  aolution  is  the  only 
solution  to  (3)  with  6 ■ 0,  note  however,  that  for  c 0 the 
solution  timy  not  be  unique. 


The  perturbation  an  tloncr.ibcd  licfL'  ia  n I’.oiuopfunl  ili'Vj(!ii!  It 
1b  linplemontbJ  in  a computer  code  by  inaklni’i  lexlcograpljlc  compariannH , 
BQQ  Dantelg  12,  Chapter  10], 

The  next  step  of  the  algorithm  is  to  "complementary  pivot"  on 
(5)  i la  Lomkn  wheroln  one  begins  with  the  initial  solution  and 
increases  0 . To  he  more  explicit  here  we  need  some  definitions. 

Due  to  the  perturbation  every  solution  of  (5)  has  (m  + J.) 
or  (m  + <,  + 1)  positive  components.  If  a solution  has  exactly 
m 1-  £ positive  components  it  is  defined  to  be  a basic  solution  and 
the  positive  variables  are  called  the  basic  variables,  by  a ray  of 
solutions  to  (S)  we  mean  a family  of  Bolutions  to  (.S)  of  form 
(a*,  t*,  X*,  p*,  0*)  + r<8,  t»  X,  p,  0)  where  r varies  over  the 
interval  [0,  +«)  and  (a,  t,  X,  p,  o)  is  nonsero. 

Stop  of  the  nlgotlthm  is  to  complementary  pivot;  such 

is  begun  In  t to  ration  1. 


Itoratlon  k • Is  Heginnlno  with  the  initial  aolutlon,  which 
is  basic,  increase  6 and  adjust  basic  vnriobles  in  order  to  retain 
a solution  to  (5);  either  wo  can  increase  6 to  infinity  in  this 
manner  or  some  basic  variable,  now  designated  the  bJ^ockijn^  variable, 
is  driven  to  eero.  In  the  first  case  a ray  of  solutions 

*’e'  ^E*  'V  generated 

where  the  positive  compononts  of  each  solution  are  0 and  the  basic 


9 


vnriablpBi  the  pivotlna  la  termlniiCo.il . In  the  aucoiid  iwum  a new 

solution  , , 11^,  G^l  tins  been  generntod  wherci  tlio  poalUve 

i c e e t t.  j 

components  are  6 and  the  basic  variables  excluding  the  blocking 
variable:  the  new  aolutlon  is  basic  and  we  move  to  Iteration  2.  H 


UcxnjtJijm 

(k+l)  ■ 2 

. 3 , > , , ! 

Asaume  that 

tho  haalc  aolutiops 

/ 0 

I B,  1 

t”, 

,0  0 

1'  • 

o”),  .... 

.k  k nk^ 

have  been  generated 

V " 

L 

C F, 

y 

1 c'  e 

i;  t f.  y 

and 

that 

oj^  la 

positive. 

\ 

There  is 

exactly  one 

1 In  {1 m+A) 

such 

that 

bo  til 

and 


0 . 


One  of  those  two  variables  was  basic  in  the  basic  solution} 

let  UH  now  doslgnatu  the  other  variable  of  the  pair,  the  drlvl na 
varigblo.  Beginning  with  the  basic  solution  X^,  ^1"', 

increase  the  driving  variable  and  adjust  the  basic  variablea  to  main- 
tain a solution  of  (5);  either  wo  can  Increase  the  driving  variable  to 
infinity  in  this  manner  or  some  basic  variable,  now  duaignntcd  tho 
blocking  variable,  Is  driven  to  zero,  In  tlie  first  case  a rny  of 
solutions  ^8^,  tp,  x|^,  u^,  + r^s,  t,  I?,  JI,  oj  to  (5)  has 

been  generated  where  the  positive  components  of  each  aolutlon  are 
tho  driving  and  basic  variables;  the  pivoting  la  terminated.  In  tits 

second  case  a new  solution  been 

I c e ' E ’ c t y 
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gt'.ncrnlod  wluirr  t:lio  poHllive  lUMiipmu'iil  m .-iro  1.110  (Irlviiij;  ,'iikI  h.isic 
vaiii.flblc’H  cxo  I iicllii!’  l ilt?  blocking,  v.-.ii:.Liib.l.o;  tlio  tuw  tio,lut:.ton  Im  Ixifilc 
and  we  move  to  Iteration  k + 2.  M 

Our  di'ftcr.lpt.ion  of  the  nlfioritlim  ia  now  oHSontially  complo.tc! . 

The  only  remaining  ni.TtLer  la  to  show  that  it  generated  the  path  of  the 
Theorem. 

We  lianteii  to  remark  that  0^  la  positive  for  k exceeding 
isero,  because  the  initial  solution  is  tlui  only  one  wJ.tli  0 equal  to 
zero  and  because  the  algorithm  cannot  cycle;  this  argueiuent  is  ampl.lCied 
upon  in  Appendices  1 and  2.  Since  there  are  only  finitely  many  baaic 
solutions  the  complementary  pivoting  must  terminate  on  a ray. 

Assuming  the  complementary  pivoting  terminates  in  iteration 
k 4-  1 , the  ray  of  solutions  to  (5)  generated  lia.i  the  form 

^e‘  'V  °c)  where  (s,  t,  A,  ii,  t))  / 0 

and  r varies  over  (0,  . 

Define  x*^  and  x by 

E 


is  n rny  of  .so.lutionH  to  (6)  wlicro  r vnfios  ovtu'  tho  iiiLovvol 

[0,  •+’'-)  . 


Of  course',  It  fo]]ows  that 


(7) 


X 8 t A n 0 


c 

0 

0 

a''' 

u''' 

0 

A 

1 

() 

0 

0 

0 

U 

0 

1 

0 

0 

s>0  t>0  A > 0 |i>0  U>() 


n*A  ^ = tni*^  “ 8^*A  '■  h*^*A*^  k 

c t;  e e i:  r.  t t 


Towards  provlnp  our  tiiaJn  thuorem  it  is  Important  to  prove  that. 

5 1h  poHitivo.  If  the  complomontary  pivoting  torin.luatLiH  in  iteration 
one,  tluMi  obviously  0 is  positive.  If  the  pivoting  t orml ootns  in 
itorntlon  k + 1 where  k exceeds  one,  then  0^  is  positive  as 
we  have  aircady  argued.  I.et  u.s  suppose  0 is  aero  and  show  that  it 

leads  to  a contradiction. 

Suppose  0 iij  aero.  Then  x is  zero,  since  that  is  the  only 
solution  to  Ax  0 and  Hx  £ 0 . Hero  a and  t are  zero.  Thore- 
■fore,  + fill  ■ 0.  Adding  ^Ax^  “ A(a  + [t]^)  to  pUxJ^  - iieoj' 

- p(h  -t  [l]^)  we  get  -pi'O^  « A(a  + [c.]^)  + p(b  + ft  1^)  . So  If  u 
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iWu  I 


c 


J 


Is  non/.oni,  t.lii.'  limit  i.pJ.f.orn  {X,m)  show  Ax  ;',  a |i;],^  with 
lix  < li  + to  be;  Inroniiib.l  C!  wlilrlt  i.s  a (innt  nid  I i-l  I on , Tlui.s  , p 

is  zero  and  A Is  nonzero.  IkMico,  AA  0 and  A (a  -I-  (i  ]^)  ™ ()  , 
it.  ioHows  tliat  (lt]..)j  zero  whtuiov'or  A,  Is  pos.ltlve.  lUil 

J J 

this  cannot  be  by  t'ho.lce  of  (c]2f  alno.e  tlio  rows  A.^  ^ for  A^ 
po.slLlvo  arc  linoarly  di'itinulem  . Our  .siippi>«.l tion  .l.omln  I'o  a r.nnt.ru- 
d.iotion  and  wo.  m.iy  a.4nnme  0 :l.s  po«.It;lvr. 

Now  tbiif  0 J.H  OHtabllahod  as  pos.lt.tve,  wo  ro, lax  our  portur- 
bation,  tliat  lu,  for  .i  “0,  .1,  k l.lio  ill)  so.liitlon 

(x^,  A^^,  0^*)  to  (fi)  bo.ooiiios  a solution  (a*,  s*,  A^,  p\  0'^) 

to  (2)  upon  .sottlnp,  i'  to  zoim;  a co(ird'l.nat:o  of  flu;  lator  solution  Is 
po.sitlvo,  only  'if  the  (‘orro-spondlnc,  solul  Ion  of  t lie  ft'i'iiiur  one  Is 
pof.tI tlvo.  Wo  have  that  (x*^,  t*^,  A*^,  ii*^,  0*^)  -I-  r(M,  s,  t,  A,  ii,  ii) 

ns  a ray  of  solutions  to  (2). 

To  prove  tlio  'I'lioorem  wo  sot  X(l)  and  0(i)  “ 0*  for 

1 “ 0,  ki  oxtond  (N.O)  to  (0,k]  by  makiiij;  It  nffliu;  on  oath 

Interval  |l,  1 + IJ  fiU'  I 0,  k-i,  and  finally,  oxtond  (X,0) 

to  [0,  +"')  by  .soitlnp,  X(l(+r)  “ + rx  and  O(lol-r)  ■"  + rti  for 

all  r > 0 . 

k k 

To  prove  tlio  Corollary  we  first,  consider  p , If  p Is  zero 

then  X Is  a stationary  point  of  (Jf.JC)  . Othorwiso  let  us  aasumo 

k k - 

U Is  not  zero  and  wo  coiuitder  the  ray  x + rx  in  JC*  Wu  have 
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X • + rx)  “ -X  • t-  rA)  - x • + I'lO 


s • (\**^  + rA)  + (t  - 0l‘)  • (|i^  4'  riO 


“ -Oo  • (ii'^  4 rii)  “Oo  • ii^  < 0 


and  thn  ri>m»1.l  Is  onliiblished  whom  x*  x'^  . 


3-  CW  >’).  << 

l,ti  t'.lu'  provloiiH  nuctlon  flu'  quantitlisfl  ()J,  b)  nnd  fi 

wuro  tiGctlcil  Lo  execute  the  alpjorltlim.  In  pafticulnr  Hott.Inj’s  tlieli' 
eelcotlon  mJp,lil;  bo  evident,  for  example,  in  fX  minl't  be  n 

known  ositlinato  of  a atntlonary  point.  Our  piirposio  here,  howovor, 
if?  to  doMcrIbu  a Koneral  way  of  K‘UK*rat In(i  x^,  (II,  b)  » and  « . 

Soiect Lnc  x^*i  Uho  Phase  T of  the  fiimp.lox  metliod  to  solve 
tile  ayatem  Ax  + la  " a,  a > 0 . Anfluinlii}(  Ax  < a is  fonnible, 
the  system  (H)  in  obtained  via  elementary  row  opera l;Jonn  on 
Ax  + IB  “ il  . 


Note  that  0,  that  iJ,  n)  is  a disjoint  union  of  y and 

6,  that  {],  . . . , ra)  is  a dlnjolnt  union  of  ti  nnd  f,,  and  that 

Y nnd  n contain  the  sarao  nuinbe.r  of  clementH.  (Henri y x^  defined 

by  x^  “ Oj  and  x^*  “0  is  a solution  to  Ax  < o . U) 


I 


I 

i ^ 


So  1 1 1 nj;  (!',  I))!  fSY.'ilo 


A X 

H* 


(9) 


’‘6  - ^ 


o • A X + ox.,  >-■ 

ti  • 6 - 


- 0 • II 


lias  a un.l((iip  solution,  namoly,  x^  , To  soci  thin  obnurvo  that  (9) 
is  or|iilva.tc*nt  to 


A X »*  a , 
n*  n 


x^  " 0 


ami,  conHciiuontly , to 


:i 


I|  .1 

:j  i 


A X “n  a , 
iiY  Y n 


x^  - 0 . 


But  from  tlio  olomontary  row  oiHirnt.iovifl  wc>  see  that  A„  a « a,  . 

n y n 1 

Thurc'Foro,  wo  can  (lof.lnti  (U,  h)  by 


-I 


rt  • 


0 1 ,,  - oA 

k)  • n • 


b - 


-o  • a 


•4  . 


S(?lc!c_tln('  h!  I-ot  [i  be  tho  set  n tJ  {m  + 1,  in  + £}  of 

n + 1 elements  whore  B is  ii.  x n . Solve 


(' 


T T 
a‘,  I)‘ 


’P  v'j  P 


0 

-Cx  - c 


(X,  tOji  1 0 


for  n hnnlc  solution.  Then  discard  from  p an  element  corro.spondliiR 
to  a Koro  (’omponoiit  of  llui  nolution  to  ('ot  B . IS 
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. kxamjm  ,h: 


ll,si»n  tlu'  u^f'orlt'.liin  rirj  Mpoci  f lc'<l  dn  Suctionii  2 nn<1  3 wc  com- 
pute H Htntionnry  point  of  (j^’,0^)  whom  ^(x)  - Cx  + {x  ! Ax  < n} 

(tnd 

/ ,1  0 0 2 

(A,  n)  - j -1  2 0 6 

^-2  -A  0 -A 


By  proJocUiiiR  ^ Into  the  (xp  x^)  plone  wo  p.ot  tho  hoc  of  FiBure  I. 


UlKin  fi])plyin(i  l’h(i,si;  I to  t.lu!  .sviUi'iii  Ax  -I-  l;s 


il  , H ■ 0,  til.  1 1 I'i. 


tlio  f^yHtc’nl 


’‘1 

"3 

”1 

"2 

^3 

I 

0 

] 

0 

.1 

1 

i ^ 

*'  1 

' 

! 0 

(1 

-2 

-4 

0 

0 I 0 

‘ 

-4 

HyHtt'in 

(R) , iiiiiiiely 

^,1 

^3 

®2 

”3 

1 

0 

n 

:i 

0 

0 

?. 

-2 

0 

0 

0 

1 

.5 

4 

.3 

.1 

0 

0 

0 

-.r-is 

1 

wlicroupoii  we  get  - (0,  1,  0) , 7 - { 2) , i3  - U , .'3) , f,  - {] , 2 ) , 
n “ { 3 } , S ■ { 3 , 4 , 5 , 6 ) , a nc* 


- — MuxwBmpujim 
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BESI  AVAILABIE  COPY 


Plmir<’  2 tllfipliiyn  HlniultanooiMtly  Hy(il:('nis  (A)  aiKl  (S) . lA)i' 

a 

? 6 I ■' 

a glvon  row  tlic.  port;iirbatl.oii  coo>  fli- l.oiil'i)  f»r  t. , r.  , i;  an' 

found,  renpool :l voly , In  columna  f),  10,  12,  1,  2,  .'3;  (IiIk  fiiul  rwiialim  | , 

It 

!io  following  plvol'.ii  on  ctio  matrix.  ' | ■ 

By  i))volinR  on  the  poKitlons  (1,  9),  (2,  10),  and  (3,  12)  of  | 

tlio  matrix  of  I'iguro  2 tlto  Inllinl  riolution  (k  " 0)  iH  dl.nplayod,  } ' 

Uliat  1h,  luivo  tlio  cminonlf.al  form  with  rcupoc. t to  the  inJllnl  | ■. 

I. 

haaic  Bolntlon.  Next,  for  tho  c.omploinuntary  pivoting  wo  continuo  i, 

) • 

pivoting  on  tho  matrix  on  poHltlona  (1,  13),  (6,  3),  (S,  11),  (2,  8),  ; 

(5,  A)  and  (3,  li)  in  order  to  oxoeuto  Iterations  k " 1 , 2 , . , . , fi  i 

wharonpoti  a ray  Is  cncoun tercel.  The  aolufloiuj  corrunpondlnn  to  each 

I . 

Itnratlon  are  glvcMi  in  Figure  .3i  In  ooiiuiinn  h tlio  values  for  i 

(x,  ii,  t,  1,  ii,  0)  aro  diaplaycd.  Aa  antic Ipatecl  from  tlio  Coiollnry  | 

the  nlgorlthm  yields  a stationary  point  for  (jfj,9c) , namely, 

X - (-1,  2.5,  -1)  . In  Plguro  A the  path  (X,0)  an  pro, lee  led  onto 
the  (Xj , x^)  space  is  shewn. 


APPENDIX  1;  UNIQUE  INITIAL  SOLUTION 

The  following  lemma  shows  that  (5),  or  equivalently,  (6), 
has  a unique  solution  with  6 « 0 for  all  sufficiently  small  positive 
e . Observe  that  it  Is  critical  that  the  constraints  are  perturbed 
less  than  Here  our  notation  has  been  simplified,  and  the  result 
Is  more  general  than  is  required  in  Section  2. 

Let  be, the  set  of  x's  in  r"  that  satisfy  Ax  _<  a + BteJj 

where  B la  m x m and  [e]2  ■ t”"*™)  ► Let  be 

the  linear  function  Cx  + c + D(e]j^  where  D is  n x n and 
[e]^  ■ e")  . Consider  the  linear  program  (y,  e) 

min  X • ^j.(y) 

s/t  Ax  ^ a + B[e]2 

Of  course,  if  x is  an  optimal  solution  to  the  program  (y,  c)  and 
if  y “ X,  then  x is  a stationary  point  of 

Lemma:  Assume  Ax  ^ a has  a unique  solution  x^  and  that  D Is 
nonsingular.  Then  there  Is  an  X so  that  x®  + Xte]2  Is  the  unique 
solution  of  the  linear  program  (y,  e)  for  any  y in  for  all 

sufficiently  small  positive  c. 


25 


» 


Proof;  Consider  the  linear  program  (x^,  e) . For  each  (x°,  e) 
there  is  a e such  that  exists,  x - A“^(a  + is  optimal, 

and  fo  (x^)  < 0.  Consequently,  since  there  are  only  finitely 

P* 

many  3's  and  D is  nonsingular,  there  is  a 6 such  that 

a - 3 and  (a"^)"^  ^ (x°)  < 0.  for 

p c 

all  sufficiently  small  positive  e.  Next  observe  that  e times  the 

diameter  of  ^ tends  to  zero  as  c tends  to  zero.  Hence,  for  all 
£ 

sufficiently  small  positive  e and  any  y in  we  have 

A^(A"^(a  + B[e]2)g)  1 and  ^^(y)  < 0.  It  then  follows  that 

for  all  sufficiently  small  positive  e that  x^  ■ '*■ 

the  unique  optimal  solution  of  the  linear  program  (y,  e)  for  any  y 

in  if  x / is  in  then  A^.x  !^(a  + B[e]2)^. 

( <?g(y))'”’  A^^  Ag.x  > ( ]&g(y))'^  A“J(a  + B[e]2),  and 


^,(y)  >.x^ 


(y).  ta 


T 


(■onvi'.ucknci-:  to  a ra\ 

Till})  iipiM'iulix  1h  i-iuiipJ  it'd  to  assist  t lie  rc-niifi'  who  iU',:(ri  s to 
pursue  more  ccit'fully  tliu  (;oiivov(>onco  proof  I'or  Ltio  c’oir.i  r.ioiitury 
pIvot.liH’,  of  Sof't'lon  Llfilnn,  tVio  ooncopts  of  [51  wo  atpou.',  briefly, 
that  t.lu.'  filp.orJ  t.lim  noiioriiles  a r;iy.  A thorouf'h  arj;ui'i''Ut  Is  not  [iisti- 
flod  hove  lio.o.niHo  It  is  lonf.,rliy  and  booiiiiue  the  1 1 ' oi'a |■.ul•l.>  la  now 
laden  with  proofs  est)tbll  alilnjt  .similar  conohuiloi,  : l.'wovt'r,  uainft 
tl\a  rein, arks  mado  belciw  Lot'othoc  with  tlio  p.onoriil  t'u'oi'y  of  [.5],  one 
has  a rigorous  and  coiiiploto  wxpositlon  of  tlio  o.(>nvor)',tMu;o  lo  a ray. 

As  in  Appoiulix  1 wo  simplify  the  notation. 

Con.'.!  id  or  the  systo.m 

Uw  + Sx  -h  dQ  « q + Q[f.1 

(10) 

w > 0 2 > 0 0 > 0 w • z " 0 

where  the  in.itrtui'.s  R,  S,  and  Q are  n x n,  Q Iss  nonslnpular , the: 

vectors  w.  z,  d,  q and  (c. ] are  n x ],  [t]  *•  (i  \ and 

the  vnrlal.i.l.eu  are  (w,  z,  0)  . Identify  this  tiyatoin  wl  tli  f.'i)  where 

w - (s,  t),  z “ (X.li),  q " (q^,  q^) , etc. 

In  Appendix  1 we  proved  that  for  each  .sufficiently  airiall 
positive  i;  the  systom  (5)  has  a unique  solution  with  0 zero,  and 
now  let  us  nssuine  that  (10)  also  haa  this  property. 


BtSl  AVMlABLt  COPN 


l.et  M 1^C:^  1;1k'  tuM.  R x /nul  let  0'f\  bo  the  ^^ub(l.l vinlon 
of  M will  MO  col  Is  1)1  nro  of  form  C/  x .■ind  ulu  lo  C'^  In  nny 

orLhnnt.  of  u'’  . l)ef:l.tiu  the*  pJccovv^ifit*  linear  map  I'  : M + K*'  by 

n 

F(y.0)  " I 1%  (y.)  + Del 

1 ^ ^ 


where  defined  by 


if  j*^  0 
If  y^  < 0 


i 

b 


First  wo  observe  that  the  syotum  (J.O)  In  equivalent  Lo  F(yi  0)  * q + [(.  ], 

and,  in  particular,  F(y,  0)  " q + [r]  has  n unique  fsoiution  for 

•"3  r\ 

small  positive  e,  or  in  other  v-’ordn,  (F  (q  + [t])  (‘1  (l<  x 0)  con- 
tains exactly  one  oleinont  whore  r"  x 0 is  the  boundary  of  M . 

For  small  positive  c,  F ^(q  + [c])  is  a 1-mnnlfold  neat  in 
(M,^t  let  W^,  be  the  route  in  F ^ (q  [e])  which  mnetn  r”  x q . 

Hie  complementary  pivot  scheme  applied  to  (10)  follows  for  all 

sufficiently  small  t.  be);innlnj’  with  the  point  i i (It”  x o)  . 
is  subdivided  by  a finite  number  of  cells  of  form  Ho  where,  o 
la  an  element  of  Since  cannot  return  to  the  lioundary  of  M, 

it  must  terminate  with  a ray. 
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